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Abstract 

We discuss low-temperature multipole states of Nd-based filled skutterudites by analyzing a multiorbital Anderson model with the 
' use of a numerical renormalization group method. In order to determine the multipole state, we take a procedure to maximize the 
multipole susceptibility matrix. Then, it is found that the dominant multipole state is characterized by the mixture of 4u magnetic 
and 5u octupole moments. The secondary state is specified by 2u octupole. When we further take into account the coupling between 
/ electrons in degenerate Y§ 7 (e u ) orbitals and dynamical Jahn- Teller phonons with Eg symmetry, quadrupole fluctuations become 
significant at low temperatures in the mixed multipole state with 4u magnetic and 5u octupole moments. Finally, we briefly discuss 
possible relevance of the present results to actual Nd-based filled skutterudite compounds. 
[ © 2008 Elsevier B.V. All rights reserved. 
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1. Introduction 

Recently, magnetism and superconductivity of rare-earth 
and actinide compounds have attracted renewed attention 
in the research field of condensed matter physics [1]. In 
particular, filled skutterudites, expressed as RT4X12 with 
rare-earth atom R, transition metal atom T, and pnictogen 
X, provide us a platform for systematic research of mag- 
netism and superconductivity of /"-electron systems with 

\ n>2 [2,3], where n denotes the number of / electrons. 
Since RT4X12 crystallizes in the cubic structure with 
high symmetry of Th point group [4] , orbital degeneracy re- 
mains in general. Due to the strong spin-orbit coupling in 

I / electrons, spin-orbital complex degrees of freedom, i.e., 
multipoles, become active in filled skutterudites. For in- 
stance, a second-order phase transition at 6.5K in PrFe4Pi2 
[5] has been considered to be due to antiferro quadrupole or- 
dering [6] . Note that a possibility of antiferro hexadecapole 
order has been also suggested in PrFe4Pi2 [7]. In NdFe4Pi2, 
a significant role of quadrupole at low temperatures has 
been suggested from the measurement of elastic constant 
[8]. A possibility of octupole ordering in SmRu4Pi2 has 
been also pointed out from the elastic constant measure- 
ment [9]. Note that the octupole scenario in SmRu4Pi2 
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has been supported by muon spin relaxation [10] and 31 P 
NMR experiments [II] . Quite recently, a possibility of anti- 
ferro hexadecapole order has been proposed to understand 
metal-insulator transition of PrRu4Pi2 [12]. 

Another characteristic issue of filled skutterudites is 
rattling, i.e., anharmonic vibrations of rare-earth atom 
around the off-center position inside the pnictogen cage. Ef- 
fects of rattling on low-temperature /-electron states have 
been recently discussed actively, in particular, with rele- 
vance to magnetically robust heavy-fermion phenomenon 
observed in SmOs4Sbi2 [13,14]. Concerning the symme- 
try of vibrations, a possibility of degenerate E g mode has 
been suggested in PrOs4Sbi2 [15]. Since there exists linear 
coupling between / electrons in degenerate r^ 7 (e u ) or- 
bitals and vibration mode with E g symmetry, the present 
author has pointed out quasi-Kondo phenomenon due to 
dynamical Jahn- Teller (JT) phonons [16,17]. 

In this paper, we focus on the case of n=3 as typical ex- 
ample to study low-temperature multipole properties and 
the effect of JT phonons on the multipole state of filled skut- 
terudites. The multiorbital Anderson model constructed 
based on a j-j coupling scheme is analyzed by a numeri- 
cal renormalization group method. Note that the multipole 
state is determined by the maximization of the multipole 
susceptibility. It is found that the primary multipole state 
is characterized by the mixture of 4u magnetic and 5u oc- 
tupole moments, while the secondary state is specified by 
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2u octupole. When we further include the coupling between 
/ electrons and JT phonons, we find that quadrupole fluc- 
tuations are significant at low temperatures in the 4u-5u 
coupled multipole state. Finally, we briefly discuss possible 
relevance of our results to actual Nd-based filled skutteru- 
dite compounds. 

2. Multiorbital Anderson Model 

The local /-electron state is described by [18,19] 



Hi, 



{^m^m'^crcr' H" ^Cm^o^m' ,(7') frna fm' 



^mi,m2,m3,m4/ miffl /m 2 iT 2 /m3iT2/m4ffi > (1) 



777,(7, 771 ' .(7 

EE 

mi~m4(Ti,(72 



where f ma is the annihilation operator for / electrons with 
spin a and angular momentum m(=— 3,- • -,3), <r=+l (—1) 
for up (down) spin, B miTn i is the crystalline electric field 
(CEF) potential for angular momentum £=3, 8 aa : is the 
Kronecker's delta, and A is the spin-orbit coupling. The 
matrix element Cm,cr,m',cr' is given by 



m,±l,m,±l 



±m/2, 



Cm±i,=Fi,m,±i = \/l2- m(m ±l)/2, 



(2) 



and zero for the other cases. The Coulomb integral 
^TTii, ma, 77i 3 , 77i 4 i s expressed by the combination of four 
Slater-Condon parameters, F°, f 2 , f 4 , and F 6 [19]. In 
this paper, we set F°=10, F 2 =5, F 4 =3, and F 6 =l in the 
unit of eV. For the Th point group, -B m , TO ' is given by three 
CEF parameters, B%, B%, and i?| [4,18]. In the traditional 
notations [20,21], they are expressed as 

Wx/15, B° = W(l - |£c|)/180, B\ = Wy/24. (3) 



where x and y specify the CEF scheme for Th point group, 
while W determines an energy scale for the CEF potential. 

The local Hamiltonian H\ oc can provide us exact infor- 
mation on local /-electron states, irrespective of the values 
of Coulomb interactions and spin-orbit coupling [18]. How- 
ever, since H\ oc includes seven orbitals, we are immediately 
faced with difficulties for further study of many-body phe- 
nomena in /-electron systems. Thus, it is natural to con- 
sider the effective model which describes well low-energy 
states of H\ oc . For the purpose, we have proposed to ex- 
ploit a j-j coupling scheme [22,23]. We set the spin-orbit 
coupling term as an unperturbed part, while the CEF po- 
tential and Coulomb interaction terms as perturbations. 
Then, we obtain the effective model of H\ oc as [23] 

-ffirfF = ^""^ B ft, v fl_.fi; + y ] ^i,^2,^3,M4//ii/ii 2 /f [i 3/f[*4; (4) 



where / M is the annihilation operator for / electron with 
angular momentum 5/2,- • -,5/2) in the j= 5/ 2 sextet. 
The modified CEF potential is expressed as 



where B^} v denotes the CEF potential for J—h/2 and B^]l 
is the correction in the order of W 2 /X. The effective inter- 
action in eq. (4) is given by 



1 Ail,/J2,/J3,M4 



where ij^i,n 2 ,H3^i is expressed by three Racah parameters, 
Eq, Ex, and E 2 , which are related to the Slater-Condon 
parameters. Explicit expressions of by using Eq, E\, 
and E 2 are shown in Ref. [22] . 

On the other hand, ^1^2,^3,^4 is the correction term 
in the order of 1/A. Details on this term have been dis- 
cussed in Ref. [23]. Here, three comments are in order, (i) 
Effects ofBg and B% are included as two-body potentials in 
/W. (ii) The lowest-order energy of /W is \W\J K /\, where 
Jh denotes the original Hund's rule interaction among / 
electrons, (iii) The parameter space in which H c ff works 
is determined by the conditions for the weak CEF, i.e., 
|W|/J H <1 and \W\ J H /A<C_E 2 - Since E 2 is the effective 
Hund's rule interaction in the j-j coupling scheme, esti- 
mated as £ 2 ~Jr749 [22], we obtain |W|/A<0.02. Thus, it 
is allowed to use H e g even for A in the order of 0.1 eV [23], 
when \W\ is set as a realistic value in the order of 10~ 4 eV 
for actual /-electron materials. 

Now we consider the hybridization between / and con- 
duction electrons. From the band-structure calculations, it 
has been revealed that the main conduction band of filled 
skutterudites is a u with xyz symmetry [24], which is hy- 
bridized with / electrons in the state with a u symmetry. 
In order to specify the /-electron state, we introduce "or- 
bital" index which distinguishes three kinds of the Kramers 
doublets, two Tgy and one . Here "a" and "b" denote the 
two r^" 7 's and "c" indicates the . 

Then, the multiorbital Anderson model is given by 

^ = E £feC L c ^ + E( yc L/c CT + h.c.) + H cS + H cph , (7) 



(Ul,/i2,(U3i(I*4 ' AH,(U2,(li3,/i4' 



(6) 



ka 



ha 



b - B (0) 4- J5 (1) 



(5) 



where is the dispersion of a u conduction electrons with 
r^" symmetry, f ia is the annihilation operator of / electrons 
on the impurity site with pseudospin a and orbital 7, cuo 
is the annihilation operator for conduction electrons with 
momentum fc and pseudo-spin a, and V is the hybridiza- 
tion between conduction and / electrons with a u symmetry. 
Throughout this paper, we set ^=0.05 eV. Note that the 
energy unit of H is half of the bandwidth of the conduction 
band, which is considered to be in the order of 1 eV, since 
the bandwidth has been typically estimated as 2.7 eV for 
PrRu4Pi2 [25]. Thus, the energy unit of H is taken as eV. 
To set the local /-electron number as n—3, we adjust the 
/-electron chemical potential. 

The last term in eq. (7) denotes the electron-phonon cou- 
pling. Here, the effect of E g rattling is included as relative 
vibration of surrounding atoms. We remark that localized 
Tgy orbitals with e u symmetry have linear coupling with 
JT phonons with E g symmetry, since the symmetric repre- 
sentation of e u xe u includes E s . Then, i/ op h is given by 
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H cph = g(Q 2 T x + Q 3 t z ) + (P| + P 3 2 )/2 

+ (uj 2 /2)(Q 2 2 + Ql) + b(Ql-2Q 2 2 Q 3 ), (8) 

where g is the electron-phonon coupling constant, Q2 and 
Q 3 are normal coordinates for (x 2 — y 2 )- and (3z 2 — r 2 )- 
type JT phonons, respectively, P 2 and P3 are correspond- 
ing canonical momenta, t x = J^MLfba + fLf^a), t z = 
- flafba), w is the frequency of local JT 
phonons, and b indicates the cubic anharmonicity. Note 
that the reduced mass of JT modes is set as unity. Here 
we introduce non-dimensional electron-phonon coupling 
constant a and the anharmonic energy (3 as a—g 2 /(2uj 3 ) 
and /3=&/(2w) 3 / 2 , respectively. 

3. Multipole Susceptibility 

In order to clarify the magnetic properties at low temper- 
atures, we usually discuss the magnetic susceptibility, but 
in more general, it is necessary to consider the susceptibil- 
ity of multipole moments such as dipole, quadrupole, and 
octupole. The multipole operator is given in the second- 
quantized form as 

A 7 = J2(X,),„flU, (9) 

where X denotes the symbol of multipole with the symme- 
try of T 7 and 7 indicates a set of indices for the irreducible 
representation. For j—5/2, we can define multipole opera- 
tors up to rank 5 in general, but we are primarily interested 
in multipole properties from the T^ 7 quartet. Thus, we con- 
sider multipole moments up to rank 3 in Oh symmetry. 

Now we show explicit forms of multipole operators 
[26,27]. As for dipole moments with T 4u symmetry, the 
operators are expressed as 

Jaux — Jxi ^4uy — Jyi Jiuz — Jzi (10) 

where J x , J y , and J z are three angular momentum oper- 
ators for j=5/2, respectively. Concerning quadrupole mo- 
ments, they are classified into r 3g and r 5g . We express the 
T3 g quadrupole operators as 

3g « = (2J 2 - J 2 - J 2 )/2, 3gv = V3(J 2 - J 2 )/2. (11) 
For the r 5g quadrupole, we have the three operators 

C>5g£ = V3JyJ Z /2, 

5gv = V3XX/2, (12) 

05g£ = VS J X Jy/2, 

where the bar denotes the operation of taking all possible 
permutations in terms of cartesian components. 

Octupole moments are classified into three types as r 2u , 
T4 U , and r 5u . Among them, r 2u octupole is written as 

T 2u = Vl5J x JyJ z /6. (13) 



3 

For the r 4u octupole, we express the operators as 



T^uX (^J X J X Jy JxJ Z )/2, 

T 4uy = (2J a 3 -Jpf-Jpf)/2, (14) 

T^uZ — (2J Z JzJ X JzJy)/2, 

while r 5u octupole operators are given by 
T 5ux = VTsiTJ 2 - Xl 2 )/^ 

T 5uy = Vl5(j y Jl-J y -Jl)/6, (15) 
T 5uz = VIE(Jz~J 2 -Jz~J 2 )/6- 

Note that we redefine the multipole moments so as to satisfy 
the orthonormal condition Tr(A 7 AT 7 /)=<5 77 ' [28]. 

In principle, the multipole susceptibility can be evalu- 
ated in the linear response theory [27], but we should note 
that the multipole moments belonging to the same sym- 
metry can be mixed. In order to determine the coefficient 
of such a mixed multipole moment, it is necessary to find 
the optimized multipole state which maximizes the suscep- 
tibility. Namely, we define the multipole operator as 

M = J2p-<X 7 , (16) 

7 

where the coefficient p 7 is determined by the eigenstate 
with the maximum eigenvalue of the susceptibility matrix, 
given by 

1 ^ e -E n /T _ e -E m /T 

X77' = 7L f _ f (n|X 7 |m)(m|X 7 /|n). (17) 

Here E n is the eigenenergy for the n-th eigenstate \ri), T 
is a temperature, and Z is the partition function given by 

In order to evaluate x 77 ' as well as an entropy Si mp and 
a specific heat Ci mp of / electrons, we resort to the numeri- 
cal renormalization group (NRG) method [29,30], in which 
momentum space is logarithmically discretized to include 
efficiently the conduction electrons near the Fermi energy 
In actual calculations, we introduce a cut-off A for the log- 
arithmic discretization of the conduction band. Due to the 
limitation of computer resources, we keep m low-energy 
states. In this paper, we set A=5 and m=3000. Note that 
the temperature T is defined as T=h~( N ~ v >/ 2 in the NRG 
calculation, where N is the number of the renormalization 
step. The phonon basis for each JT mode is truncated at a 
finite number A^ p h, which is set as A 7 p h=20 in this paper. 

4. Results 

First let us discuss the CEF states on the basis of H c g . In 
Figs. 1(a) and 1(b), we show the CEF energy levels for n=2 
and 3, respectively. Here we set A=0.1 cV, W=— 6xlCP 4 
eV, and y=0.3. Since the effects of B§ and P| are included 
in H e g as two-body potentials, H e g can reproduce well the 
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Fig. 1. CEF energy levels of H e f[ for (a) n=2 and (b) n=3. We set 
A=0.1cV, VK=-6xlO" 4 eV, and y=0.3. 

CEF energy levels of the local Hamiltonian H in the real- 
istic intermediate coupling region with A/Jh in the order 
of 0.1. As for details, see Ref. [23]. 

In order to determine the value of x for Nd-based filled 
skutterudites, here we recall that in PrOs4Sbi2, the ground 
state is singlet and the excited state is I4 2 ' triplet with 
small excitation energy as large as 10 K [31,32,33,34,35]. 
Such a situation is well reproduced by choosing x—OA for 
n=2 in Fig. 1(a). Now we change rare-earth ion from Pr 3+ 
(n—2) to Nd 3+ (n=3). In principle, it is not necessary 
to modify the CEF parameters even if rare-earth ion is 
changed, since the CEF potential is given by the sum of 
electrostatic potentials from ligand anions. Note, however, 
that the CEF potentials may be changed due to the sub- 
stitution of T and/or X in RT4X12. 

When we set x—0A for n=3 in Fig. 1(b), it is observed 
that the ground state for n—3 at x—OA is r^~ 7 quartet and 
the first excited state is T§ doublet with the excitation en- 
ergy of 0.02 eV. Experiments for NdOs4Sbi2 have suggested 
Tg 7 ground and excited states with the excitation en- 
ergy of 220K [36]. It should be remarked that the theoreti- 
cal CEF energy levels for n=3 agree well with experimental 
results for NdOs4Sbi2, by using the CEF parameters de- 
duced from the CEF energy levels for PrOs4Sbi2- We note 
that in NdFe4Pi2, both the ground and first excited states 
have been found to be Tg 7 quartet with the excitation en- 
ergy of 222 K [8]. However, in any case, the ground state 
quartet is well separated from the first excited state both 
for NdOs4Sbi2 and NdFe4Pi2- In such a situation, it is con- 
sidered that low-temperature multipole properties are not 
sensitive to the first excited state. Thus, in the following 
discussion, we fix x—OA. 

For the purpose to understand the CEF energy lev- 
els of NdFe4Pi2, it is necessary to consider first those of 
PrFe4Pi2- Here we note that the hybridization effect has 
been considered to play an important role to understand 
the difference in the CEF energy states among Pr-based 
filled skutterudites [37,38,39]. In order to determine the 
CEF energy levels of NdFe4Pi2, it is also important to in- 
clude the effect of hybridization for the case of n=3. Such 
calculations can be done, in principle, by using the effec- 
tive model H ff on the basis of the j-j coupling scheme. It 
is one of future problems. 

Now we proceed to the NRG results of H . First let us 
consider the case without the coupling between / electrons 
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Fig. 2. (a) T\ and (b) Si mp and Ci mp vs. temperature for n=3 
without JT phonons. 

and JT phonons. In Fig. 2(a), we show the multipole sus- 
ceptibility x- The dominant multipole moment in the low- 
temperature region is the mixture of 4u magnetic and 5u 
octupole, given by 



M a 



(18) 



where we find that p a =0.989, <? Q =-0.0258, and r Q =-0.146 
for a—x, y, and z. The mixture of 4u and 5u moments is 
characteristic of Th symmetry. In fact, when we calculate 
the multipole susceptibility for Oh symmetry (y=0) using 
the same parameters except for y, we actually find that 
r a —0. It is one of important features of filled skutterudites 
with Th symmetry that 4u magnetic moment is accom- 
panied with 5u octupole. Note that the magnitude of r a 
depends on parameters. The secondary multipole state is 
given by 2u octupole. We also find another mixture of 4u 
magnetic and 5u octupole moments, with reduced magni- 
tude of susceptibility. In Fig. 2(b), we show entropy and 
specific heat. At low temperatures, there remains an en- 
tropy of log 4, originating from localized r^ 7 quartet, since 
we consider the hybridization between a u conduction band 
and T§ state. In actuality, there should exist a finite hy- 
bridization between e u conduction bands and r^~ 7 states, 
even if the value is not large compared with that between 
a u conduction and electrons. Thus, the entropy of log 4 
should be eventually released. 

Next we include the effect of dynamical JT phonons, but 
before proceeding to the numerical results, let us consider 
intuitively what happens. In the j-j coupling picture, we 
accommodate three electrons into the one-electron levels 
with and r^" 7 . Note that is lower than r^ 7 , since 
the ground state of n=2 is singlet, which is mainly 
composed of doubly occupied . When we accommodate 
one more electron, it should be put into r^~ 7 . Thus, the 
Tg 7 quartet ground state is obtained. Intuitively, the 4-fold 
degeneracy is understood by the combination of spin and 
orbital degrees of freedom. 

Here we consider the JT potential in the adiabatic ap- 
proximation. Note that in actuality, the potential is not 
static, but it dynamically changes to follow the electron 
motion. For /3=0, the potential is continuously degenerate 
along the circle of the bottom of the Mexican- hat potential. 
Thus, we obtain double degeneracy in the vibronic state 
concerning the rotational JT modes along clockwise and 
anti-clockwise directions. When a temperature becomes 
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Fig. 3. (a) T X , (b) S imp and C imp , and (c) (Q z ) and ^/(Qf) (j=2 
and 3) vs. temperature for 0=0. (d) T\, (e) Si mp and Ci mp , and (f) 
(Qi) and yj (Q?) (i=2 and 3) vs. temperature for /3=— 0.002. 



lower than a characteristic energy T* , which is related to 
a time scale to turn the direction of rotational JT modes, 
the entropy of log 2 should be eventually released, leading 
to Kondo-like behavior [16], since the specific rotational di- 
rection disappears due to the average over the long enough 
time. 

In Figs. 3(a) and 3(b), we show multipole susceptibil- 
ities, entropy, and specific heat for w=0.1, a=0.5, and 
/3=0.0. Note that we suppress the cubic anharmonicity. 
We find that the results do not seem to be qualitatively 
changed from Figs. 2(a) and 2(b), in spite of the JT active 
situation. For the dominant mixed multipole moment, we 
find p a =0.902, g o =-0.412, and r a =-0.127 for a=x, y, 
and z in eq. (18). In Fig. 3(c), we show the temperature 
dependence of average displacements. At relatively high 
temperature, we find finite values of \J (Q 2 ) and \/{Qf), 
while {Q 2 ) — (Q3)—0. Namely, JT vibrations occur around 
the origin without displacements. As the temperature is 
decreased, such vibrations are also suppressed and the 
situation in the low-temperature region is well described 
by the electronic model, leading to the similarity between 
Figs. 2(a) and 3(a) at low temperatures. 

However, when we include the effect of anharmonicity, 
three potential minima appear in the bottom of the JT po- 
tential in the adiabatic approximation. Since the rotational 
mode should be changed to the quantum tunneling among 



three potential minima at low temperatures, the frequency 
is effectively reduced in the factor of e^ SE ^ , where SE is 
the potential barrier. Then, we expect to observe the quasi- 
Kondo behavior even in the present temperature range, 
when we include the effect of cubic anharmonicity. 

In Figs. 3(d) and 3(e), we show multipole susceptibilities, 
entropy, and specific heat for u=0.1, a=0.5, and /3=— 0.002. 
For T>10~ 4 , susceptibilities for both 3g quadrupole mo- 
ments are significant, suggesting that quadrupole fluctua- 
tions are dominant in this temperature region. However, 
when the temperature is decreased, \3gv is suppressed, 
while X3g« remains at low temperatures. Instead, the mixed 
multipole with 4u magnetic and 5u octupole moments be- 
comes dominant. Note that \ for M z is slightly larger than 
those for M x and M y . 

Around at T=10~ 3 , we find a peak in the specific heat, 
since an entropy of log 2 is released. As mentioned above, 
this is considered to be quasi-Kondo behavior, originating 
from the suppression of the rotational mode of dynamical 
JT phonons [16]. In this case, the entropy of log 2 concern- 
ing orbital degree of freedom coupled with JT phonons is 
released, while there still remains spin degree of freedom 
in the localized r^" 7 quartet. In fact, at low temperatures, 
magnetic susceptibility becomes dominant. 

In Fig. 3(f), we show the temperature dependence of av- 
erage displacements. For T>10~ 5 , we find \/ (QfJ^O an( i 
V (QD^Oj suggesting that both Q 2 and Qz modes are ac- 
tive. This is consistent with the finite values of suscepti- 
bilities for Ozgu and 03 g „. Note that the Q3-type displace- 
ment is considered to occur, since (Q 2 )—0 and {Qz)^Q. In 
the low-temperature region, we find v /(Q|)=(Q 2 )=0, while 
V (Qi) — I (Q3) I^Oj indicating that only Q3-type JT vibra- 
tion is active with finite displacement. This is also consis- 
tent with the result that X3gu remains at low temperatures, 
since the vibration mode is fixed as Q3-type after the quasi- 
Kondo phenomenon occurs. 

5. Discussion and Summary 

In this paper, we have clarified that the magnetic state 
with active quadrupole fluctuations appears in Nd-based 
filled skutterudites, when we consider the effect of dy- 
namical JT phonons. In fact, the existence of degenerate 
quadrupole moments has been suggested from the experi- 
ment of elastic constant [8]. In Fig. 3(d), in the tempera- 
ture region of T>10~ 4 , we have observed that both Oz gu 
and Os gv become active, although they are not exactly 
degenerate due to the effect of JT phonons. However, the 
idea of the magnetic state with active quadrupole fluctu- 
ations seems to be consistent with actual Nd-based filled 
skutterudites. Note that in the present NRG calculations, 
we cannot conclude the nature of intersite magnetic in- 
teraction, ferromagnetic or antiferromagnetic, although 
Nd-based filled skutterudites are ferromagnets. 

In Nd-based filled skutterudites such as NdFe4Pi2 [40] 
and NdRu4Sbi2 [41], peculiar behavior of a resistance min- 
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imum in the temperature region higher than a Curie tem- 
perature Xc has been pointed out. Quite recently, Np- 
based filled skutterudite NpFe4Pi2 has been synthesized 
[42]. Since actinide ion is considered to take a tetravalent 
state in the filled skutterudite structure, NpFe4Pi2 is also 
classified into the case of n—3, except for the difference be- 
tween 4/ and 5/ electrons states. In fact, NpFe4Pi2 is also 
a ferromagnet with Tc=23 K and a similar resistance min- 
imum has been observed above Tq [42] . 

In the present paper, we have observed the quasi-Kondo 
behavior for the case of n=3. When the temperature is de- 
creased, an entropy log 2 originating from the double de- 
generacy of the vibronic state is released. In other word, 
this may be quadrupole Kondo phenomenon, since quadru- 
ple (orbital) degrees of freedom are tightly coupled with 
JT phonons, as understood from Figs. 3(d) and 3(f). It 
seems to be premature to conclude the mechanism only 
from the present numerical results, but the quasi-Kondo 
behavior due to dynamical JT phonons coupled with or- 
bital (quadrupole) degrees of freedom may explain quali- 
tatively the resistance minimum phenomenon in Nd-based 
filled skutterudites. Further investigations are required. 

As mentioned in the introduction, concerning the mecha- 
nism of magnetically robust heavy-fermion phenomena ob- 
served in SmOs4Sbi2 [13], a potential role of phonons has 
been pointed out from the viewpoint of the Kondo effect 
with non-magnetic origin [14]. In this context, the quasi- 
Kondo behavior due to the dynamical JT phonons may 
be a possible candidate to understand magnetically robust 
heavy-fermion phenomenon. In fact, we have found the 
quasi-Kondo behavior also for the case of n=5, but the de- 
tails of the results on Sm-bascd filled skutterudites will be 
discussed elsewhere [43]. Here we emphasize the common 
feature between Nd- and Sm-based filled skutterudites with 
the same r^~ 7 quartet ground states. From this viewpoint, 
it may be interesting to design the experiment to detect the 
effect of rattling in Nd-based filled skutterudites. 

In summary, we have discussed the multipole state for 
n=3 by analyzing the multipole Anderson model with the 
use of the NRG method. When we do not consider the cou- 
pling between JT phonons and / electrons in Tgy quar- 
tet, we have found that the dominant multipole moment 
is the mixture of 4u magnetic and 5u octupole. The sec- 
ondary multipole state is 2u octupole. When the coupling 
with JT phonons is switched and the cubic anharmonic- 
ity is included, the magnetic ground state includes signifi- 
cant quadrupole fluctuations and we have found the quasi- 
Kondo behavior due to the entropy release concerning the 
rotational JT mode. 
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